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Abstract. The 3 He and 4 He longitudinal and transverse response 
functions are determined from an analysis of the world data on 
quasi-elastic inclusive electron scattering. The corresponding Eu- 
clidean response functions are derived and compared to those cal- 
culated with Green's function Monte Carlo methods, using realis- 
tic interactions and currents. Large contributions associated with 
two-body currents are found, particularly in the 4 He transverse re- 
sponse, in agreement with data. The contributions of two-body 
charge and current operators in the 3 He, 4 He, and 6 Li response func- 
tions are also studied via sum-rule techniques. A semi-quantitative 
explanation for the observed systematics in the excess of transverse 
quasi-elastic strength, as function of mass number and momentum 
transfer, is provided. Finally, a number of model studies with sim- 
plified interactions, currents, and wave functions is carried out to 
elucidate the role played, in the full calculation, by tensor interac- 
tions and correlations. 



PACS: 21.45. +v, 25.30. Fj 



1 Introduction 



Over the past thirty years or so, much effort has gone into trying to understand quan- 
titatively the roles that short-range and tensor correlations, and two-body components 
of the nuclear electromagnetic current play in the quasi-elastic response of nuclei at in- 
termediate momentum transfers. Yet, despite the considerable attention that has been 
devoted to this topic, many open questions remain. Complications arise, in particular, 
as a consequence of the need of (and technical difficulties associated with) providing an 
accurate description of the initial bound- and the final scattering-state wave functions, 
based on realistic Hamiltonians. 

In part, the slow progress is also due to the confusing experimental picture, particularly 
for medium- and heavy-weight nuclei, that for some time obfuscated the interpretation 
of the data. Early data [|IJ had shown that, in comparison to an impulse- approximation 
(IA) calculation using a simple (Fermi gas) model, the inclusive cross section showed an 
excess of transverse strength, mainly in the region of the "dip" between the quasi-elastic 
and the A peaks. This excess was attributed to two-body currents and 7r-production, but 
not quantitatively understood. 

The longitudinal and transverse response functions, obtained during the eighties from a 
Rosenbluth separation of experimental cross sections, seemed to indicate that, in addition, 
there was a gross (up to 40%) lack of longitudinal strength in the main quasi-elastic peak, 
and a correspondingly too low Coulomb sum rule |2|, §]. While this state of affairs is yet 
to be resolved satisfactorily, particularly for very heavy nuclei like lead, where Coulomb 
corrections are difficult j|, ||, ^, 0] there are nevertheless clear indications for A < 56 from 
the work of Jourdan , who carefully analyzed the world data on quasi-elastic scattering 
including all the known corrections, that there is no missing strength in the longitudinal 
response of medium-A nuclei (for very heavy nuclei, there are still controversial issues on 
the Coulomb corrections, see Refs. || [10], [11]]). This apparent lack of longitudinal strength 
has absorbed much of the theoretical effort of the past two decades. 

On the other hand, the excess of transverse strength — presumably due to two-body 
currents — observed in the quasi-elastic region appears to be a genuine problem. In this 
respect, the experimental situation has been put in sharp focus by the work on super- 



scaling by Donnelly and Sick [|T^] which allowed to systematically compare the longitudinal 
and transverse response functions. This work showed in the most clear way that the 
transverse strength for nuclei with mass number A=12, . . ., 56 exceeds the longitudinal 
one already in the main quasi-elastic peak by 20-40%, in addition to the excess of strength 
occurring in the "dip" between the quasi-elastic and A-peaks. This excess of strength in 
the region of the quasi-elastic peak is the main subject of this paper. The region of the dip, 
which has attracted the attention in the past, will be largely ignored as the understanding 
of this region is clouded by issues relating to pion production and the A-tail. 

Theoretical calculations of two-body contributions in the region of the quasi-elastic 
peak have been performed by many groups [T3|]- [275| using different approaches. Some of 
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these calculations find appreciable contributions, 20-40% of the transverse response, due 
to the dominant two-body terms (pion contact and in-flight, and A-excitation diagrams), 
other calculations find small, <10%, effects. It is not always clear why calculations with 
similar starting assumptions give very different results. 

In general, calculations based on an independent-particle initial state (shell model, 
Fermi gas model, possibly with RPA correlations added) give very small two-body contri- 



butions in the quasi-elastic peak jr3|-[2"0"|, |2"E| |: the pion and A terms tend to cancel 
to produce a small overall effect. The origin of bigger (20-50%) contributions to the 
transverse response as obtained in Refs. E2, 1251] is not entirely understood: the different 



treatment of the A in matter in Ref. |22j may be partly responsible p7 |. 

The model study of Leidemann and Orlandini |T7J] , in which the nuclear response was 
expressed in terms of the response of deuteron-like pairs of nuclear density, first pointed 
out that it is important to account in the initial state for the tensor correlations between np 
pairs. Only when these (rather short-range) tensor correlations were included would the 
two-body terms give appreciable contributions to the quasi-elastic response. This insight 
was quantitatively confirmed by Fabrocini |2JJ, who calculated the transverse response 
of infinite nuclear matter using correlated basis function theory including one-particle- 
one-hole intermediate states. This calculation is based on a realistic nucleon-nucleon 
(N-N) interaction and two-body terms derived consistently from the N-N interaction, 
and accounts for the interactions in both the initial and final states. It also found that 
substantial two-body contributions in the quasi-elastic peak are obtained only if the tensor 
correlations, predominantly induced by pion exchange, are retained. 



The calculation of Carlson and Schiavilla [21] was performed for 4 He using Green's 
function Monte Carlo (GFMC) techniques, a realistic (Argonne v' 8 ) N-N interaction and 
again consistently constructed two-body terms. The inelastic response could be accurately 
calculated in terms of the Euclidean response (an integral over the response function, see 
below). These "exact" calculations found that the charge-exchange character of the N- 
N interaction leads to shifts of both the longitudinal and transverse strength to higher 
excitation energies, thus producing a quenching of the response in the region of the quasi- 
elastic peak. This mechanism, however, is more than offset in the transverse channel by 
two-body currents, in particular those associated with pion exchange (required by gauge 
invariance), and hence the response is enhanced over the entire quasi-elastic spectrum. 
This enhancement was found to be substantial, and in agreement with that observed 



experimentally. The study of Ref. [f2"T]] , while providing a qualitative understanding of the 
4 He quasi-elastic response, did not identify quantitatively, however, those aspects of the 
calculation responsible for the successful prediction. 

In the present paper we study the longitudinal and transverse response functions of 
light nuclei, 3 He and 4 He, using GFMC theory. Accurate data for these responses in 
the region of the quasi-elastic peak are determined via an analysis of the world data. 
A simultaneous study of 3 He and 4 He is particularly interesting as the predicted two- 
body contributions in the transverse channel increase very rapidly between A=3 and 



2 



A=4, a feature which can give us a further handle for the understanding of two-body 
effects. The study of 4 He, including the higher momentum transfers now available, is 
especially promising, since the available data || seem to indicate that the relative excess 
of transverse strength in the quasi-elastic peak is largest for this nucleus. In order to 
include heavier nuclei and hence examine the evolution with mass number of this excess 
transverse strength, we also study via sum-rules the two-body contributions for p-shell 
nuclei, for which variational Monte Carlo (VMC) wave functions are available. 

The layout of this paper is as follows. In Sec. we discuss the determination of the 
longitudinal and transverse response functions starting from the world data on inclusive 
electron scattering, while in Sec. |3| we perform a scaling analysis in order to investigate 
the global properties of the experimental response. In Sec. |] we describe the theory 
of the Euclidean responses which link the nuclear ground-state properties to integral 
properties of the electromagnetic response as measured in (e,e'), and in Sec || present 
the model adopted for the nuclear electromagnetic current. Before comparing theory 
with experiment (Sec. ^), we carry out in Sec. |^ a model study of the relation between 
the inclusive cross section and the Euclidean response in order to better understand the 
characteristics of the latter. An extension of the study to heavier nuclei via sum rules 
is given in Sec. [|, in which the dependence of the excess transverse strength upon mass 
number is examined. In Sec. |9] we further analyze the calculated results by introducing 
various simplifications, so as to identify the most important aspects of the calculations. 
Finally, in Sec. |1(] we summarize our conclusions. 



2 Experimental response functions 

In order to determine the longitudinal (L) and transverse (T) responses, we have analyzed 
the (e,e') world data on 3 He and 4 He. A determination of the response functions in 
inclusive quasi-elastic scattering from the world cross section data has many advantages 
over the traditional approach of using data from a single experiment only. Particularly 
for medium- A nuclei the limitations of the traditional approach was partly responsible for 
the misleading conclusions mentioned in the introduction and discussed in ||. 

For the extraction of the response functions, the difference of cross sections at high- 
energy /forward- angle and at low-energy/backward-angle is used. Kinematics dependent 
systematic errors do not cancel in this difference even for measurements performed at a 
single facility, except perhaps for the errors in the overall normalization of the cross sec- 
tions. The dominant systematic errors, i.e. uncertainties in the spectrometer acceptance, 
detector efficiencies, background contributions, re-scattering, and radiative corrections are 
strongly dependent on the specific kinematics. 

To improve the determination of the response functions the difference of the L- and 
T-contributions to the cross sections has to be maximized by including data over the 
largest possible angular range. This can only be achieved by including all available world 
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cross section data. For 3 He and 4 He the use of the world data not only expands the range 
of available data in scattering angle but also increases the range of momentum transfer q 
where a separation can be done, thus leading to new information on the response functions. 



At low q the extensive sets of data for 3 He [p8[ ^9|] and for 4 He [|30| , |3~l| with good 
angular coverage have been used in the present analysis. For 3 He at high q the data by 



Marchand et al. and by Dow et al. |29| , which both cover the angular region from 90° 



to 144°, are combined with complementary cross sections by Day et al. |32| which provide 
high-energy/forward-angle data with energies up to 7.2 GeV at scattering angles of 8°. 
Similarly for 4 He the data by Zghiche et al. [30|, which cover the angular region from 75° 



to 145°, and by Dytman et al. which contribute data at 60°, are complemented with 
the forward angle cross sections by Rock et al. |33fl , Day et al. [p4| , Sealock et al. fl35| , 
and Meziani et al. [[!(J covering the angular range from 8° to 37°. 

In contrast to the analysis performed for medium- A nuclei ||, Coulomb distortions 
play a negligible role for 3 He and 4 He and no corrections need to be applied. The fol- 
lowing expression, valid in the plane wave Born approximation (PWBA), is used for the 
Rosenbluth or L/T-separation: 

E <^ £ > = ^ ^ (if = + \ (iff R ^ • (1) 

where the longitudinal virtual photon polarization e is defined as 

and varies between to 1 as the electron scattering angle $ ranges from 180 to degrees. 
Here, d 2 a/dQdu) are the experimental cross sections, uj, q and Q are the energy transfer, 
3- and 4-momentum transfers, respectively, and cr Mott is the Mott cross section. The 
structure of Eq. (|I|) shows that measurements of the cross section at fixed uj and q but 
different e allow for a separation of the two response functions Rz(q,u) and Rr(q,ui). 

In practice, the experimental spectra of the various experiments were measured for a 
given incident energy and scattering angle as a function of the energy loss of the scattered 
electron, varied by changing the magnetic field of the spectrometer. To determine the cross 
section at given values of q and u>, the data have to be interpolated. This traditionally 
was done by dividing out cr Mott from the measured cross sections and interpolating the 
responses along u/E. 

In an analysis of the world data, where the various experiments were not planned for 
an ideal coverage of the g-cu-plane, the usual scheme is unreliable due to occasionally large 
spacings between various spectra. Thus, in the present analysis an improved scheme is 
employed by first dividing out an appropriate sum of elementary electron-nucleon cross 
sections, i.e. a ep for the proton and a en for the neutron, and removing kinematical 
dependencies. Essentially what is calculated from the data is the scaling function F(y, q) 
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Figure 1: Longitudinal (□) and 
transverse (O) response functions 
of 3 He at momentum transfers 
of 300, 400, 500, 600 and 700 
MeV/c. Indicated with a + are 
the upper integration limits used 
for the Euclidean response (sec- 
tion 0). 
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du 



dudVL Za ep (q) + Na en (q) dy 
The scaling variable y is fixed by energy and momentum conservation via 



y = —q + Vuj 2 + 2uj 



m 



(3) 



(4) 



neglecting small contributions from the binding energy, the perpendicular component of 
the nucleon momentum, and the recoil energy of the residual nucleus. In the next phase 
of the analysis, the extracted F(y, q) are then used to determine F(y, q Q ) at the desired 
value q Q by interpolating F(y, q) along lines of constant y. 
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Figure 2: Longitudinal (□) and 
transverse (O) response functions 
of 4 He at momentum transfers 
of 300, 400, 500, 600 and 700 
MeV/c. Indicated with a + are 
the upper integration limits used 
for the Euclidean response (sec- 
tion 0). 



cj (MeV) 

For y < 0, F(y, q) is known to be nearly independent of q over a large range. This 
makes the present interpolation as reliable as the conventional interpolation scheme even 
if the data are separated by large values of q. For y > 0, the dependence of F(y, q) on q 
is relatively more severe, since inelastic processes contribute to the cross section. Thus, 
the q- value of the L/T-separation has been chosen to minimize the correction due to the 
interpolation at large to. 

The interpolation procedure and the separation has been tested with the data of 
the Saclay experiments [[H], |28| alone. Provided the same interpolation scheme is used, 
the published values of Ri(q,uj) and Rr{q,uj) are reproduced exactly. The improved 
interpolation scheme, using y-scaling, gives results which also are identical within the 
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statistical errors. 

With the interpolated cross section data, the response functions are extracted for 
g=300-700 MeV/c in steps of 100 MeV/c for both nuclei. The combined world data cover 
almost the full e-range, with typical values ranging from 0.05 to 0.95 for most g-sets. At 
high q this has to be compared with the results of the low energy data alone which only 
cover the region from 0.05 to 0.55. In addition, with a global analysis it is possible to 
determine for the first time the response functions at g=700 MeV/c. 

If the interpolated responses of data are plotted as a function of e, a linear dependence 
is expected. In contrast to the analysis of medium- A nuclei, in which important deviations 
were observed for high q and lu, no significant deviations were observed in the present 
analysis once the quoted systematic errors of the individual data sets are included. 

The longitudinal and transverse response functions resulting from this analysis of the 
world data are shown in Figs. |T] and 



3 Scaling analysis 

In order to show the excess strength of Rp{q,uj), in this section we study the scaling 
properties of the present response functions. Barbara et al. |37| have discussed the close 
connection between the Coulomb sum rule and the notion of ^/-scaling. More recently 
the notion of ^'-scaling was introduced by Alberico while studying the properties of 
the relativistic Fermi gas model. The application of this notion to finite nuclear systems, 
requiring the inclusion of binding effects, has been discussed by Cenni et al. [[39(1 . Guided 
by these results, Sick and Donnelly have applied ^'-scaling to a large body of inclusive 



scattering data ||12|| . Scaling in ip' has the merit to allow the study of the scaling properties 
for a combined set of different nuclei. The only relevant scale parameter in the quasi-free 
scattering regime is the Fermi momentum of the nucleus which is taken into account in 
the definition of the dimensionless scaling variable ip' (approximately given by y/kp), and 
the scaling function fijp'). 

As discussed in [|12[], -^'-scaling can also be studied for separated response functions. 
The dimensionless scaling functions fp.r are defined in [|I^| as 



fi,T = k F ^r^- , (5) 



with the factors Gl,t given in ||12|| . For the relativistic Fermi gas model and in IA, the 
universal relation 

h = f T = f (6) 

is predicted. Neglecting powers higher than two in r\p = kp/m, a relation between fp and 
the Coulomb sum rule is obtained as 

d^f RFG W = l + ^V F + --- ■ (7) 



7 



Eh 



1*- 



1.00 



0.75 



0.50 



0.25 



+ q=300 
x q=400 
o q=500 
+ q=600 
"K q=700 



1 f* 



0.00 



on 



-2 



-1 



1.0 r 



0.8 



0.6 



0.4 



0.2 



0.0 



T 



o 



+ q=300 



q=400 
q=500 
q=600 
q=700 




Xc <t> * * 
+ * 



Figure 3: The scaling functions fi and /r are shown for all q values on 
the top for 3 He and on the bottom for 4 He. The upper bands of points 
correspond to fx, the lower bands to f^. 

In Fig. |]we compare the scaling functions fii^') and /tC0') obtained for all response 
functions extracted from the global analysis of the 3 He and 4 He data. Within the error bars 
of the separated data, the longitudinal response functions scale to a universal curve over 
the entire quasi-elastic peak. Scaling of Rl(q,uj) is expected and provides a consistency 
check for the Coulomb sum rule. The results for Rx(q,uj) confirm that the basic problem 
in quasi-elastic electron-nucleus scattering is the excess strength in the transverse response. 
This excess is much larger for 4 He than for 3 He. Scaling is also observed for Rr{q,u) at 
negative values of ip', thus suggesting that processes other than quasi-free knock out can 
also lead to scaling. 
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The excess of transverse strength is particularly large for 4 He. It exceeds the longitudi- 
nal strength at all momentum transfers, and does not seem to be limited to the "dip" region, 
but affects the whole quasi-elastic peak region, extending below the 7r-production thresh- 
old. The transverse strength in the dip, which increases with increasing q, is related to 
the growing overlap between the high-energy side of the quasi-elastic peak and the tail of 
the A-peak. 

In order to study the A-dependence of this excess, we can look at the longitudinal 
and transverse responses integrated over ip' — those for 12 C, 40 Ca, and 56 Fe have been 
determined in Ref. | JL2[| . We have integrated these responses over the region of tp' that 
essentially covers the quasi-elastic peak (^'l < 1.2). When limiting the integration range 
to \if>'\ < 0.5 much of the contribution from the tail of the A is eliminated, at least for 
the light nuclei. The ratio of transverse to longitudinal integrated strength is shown in 
Fig. | 

Figure f| makes it clear that: i) the excess of transverse strength rises very rapidly 
between 3 He and 4 He, and is indeed largest for 4 He; ii) it is already large at the lowest q, 
the increase at the larger q for the heavier nuclei is mainly due to the fact that the tail of 
the A peak contributes appreciably despite the restricted range of integration in ip'. 
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4 Calculations of Euclidean response 



Since we are primarily concerned with the overall strength of the longitudinal and trans- 
verse response, we consider the Euclidean response functions, defined as [5JJ 



~ roc 

EtA^t) = / exp[-(cu - E )t] R t Aq,u) , (8) 

where the RtAQi 00 ) are the standard responses, E is the ground-state energy of the 
nucleus, and u t h is the threshold for the response of the system excluding the elastic 
contribution. The longitudinal and transverse Euclidean response functions represent 
weighted sums of the corresponding Ri(q,uj) and Rx^q.uj): at r=0 they correspond to 
the Coulomb and transverse sum rules, respectively, while their derivatives with respect 
to r evaluated at r=0 correspond to the energy-weighted sum rules. Larger values of r 
correspond to integrals over progressively lower-energy regions of the response. 

In a non-relativistic picture, the Et,l can be simply obtained from: 

E L (q,r) = (0|pt( q )exp[-(F-E )r]p(q)|0)-exp^-^|(0(q)|p(q)|0)| 2 , (9) 

Mq,T) = (0|j^q)exp[-(iJ- J E; )r]j r (q)|0)-exp^-^-^ |(0|(q)|j T (q)|0)| 2 (10) 

where the elastic contributions have been explicitly subtracted, |0(q)) represents the 
ground state recoiling with momentum q, and sums over spin projections are understood. 

In this paper we present results for the scaled Euclidean responses 

EM = ^/ T /£™ )] E LAg ,r), (11) 

[G E ,p{Q )2 )r 

where Q 2 is the squared four-momentum transfer evaluated at the quasi-elastic peak. 
This removes the trivial energy dependence obtained from scattering off an isolated (non- 
relativistic) nucleon, and the q dependence associated with the nucleon form factors. The 
longitudinal response E L (q, r) is unity for an isolated proton, and the transverse response 
Ex{q, t) is simply the square of its magnetic moment. 

The chief advantage of formulating the Euclidean response is that it can be calculated 
exactly using Green's function or path integral Monte Carlo techniques, including both 
final state interactions and two-nucleon currents. While the present calculations consider 
only A up to 4, they can be very simply extended to mass up to A=10 in direct analogy 
with ground-state calculations |4T| . In the future it may be possible to use the Auxiliary- 



Field Diffusion Monte Carlo technique developed by Schmidt and Fantoni [42] to calculate 
the response for much heavier systems. 

Other techniques have also been used to calculate the response in few-nucleon systems, 
including Faddeev methods [53, [PJ and Lorentz integral transform techniques [}55| , p) 



10 



Faddeev methods sum explicitly over the final states in the system, and hence are directly 
applicable to inclusive and exclusive responses for all possible final states. This essentially 
corresponds to a complete real-time calculation of the propagation of the system. Lorentz 
integral techniques introduce a small imaginary time component in the propagation of 
the response, directly summing over a limited region of to. For systems in which a precise 
calculation is possible, the full response can be calculated ||46|. The Euclidean response 



is the fully imaginary-time response, and hence is a more integrated quantity. While 
detailed dynamical information is more limited, it is possible to perform calculations in 
much heavier systems [47]]. 

The ground-state wave functions used in this study are obtained with variational 
Monte Carlo. They are of the general form [41| : 



i*r>= n [i-u (m 

i<j<k 



[l- £ Utofak) 



F 



(12) 



where for 3- and 4-nucleon systems |<&) is simply an anti-symmetrized product of spins and 
isospins. The central three-nucleon correlation Uo(ijk) is a scaled version of the repulsive 
central component of the Urbana-IX (UIX) three-nucleon interaction. The magnitude of 
the correlation and its range are scaled via variational parameters. The pair correlations 
Fij depend upon the pair separation and the spins and isospins of the pair: 



'.J 



r 



a 



+ u\r)Sy 



+ u aT (r)a i 



Tj + U 



tT 



(13) 



The correlation U^^ij] k) is similarly scaled from the anti-commutator part of the two- 
pion exchange three-nucleon interaction. The anti-commutator depends upon the spins 
and isospins of only the two nucleons i and j, but the spatial positions of all three. Simi- 
larly the magnitude of the spin-isospin dependent correlations u for pair ij are quenched 
by the presence of other nucleons. Both the two-nucleon correlation F^ and the Un; cor- 
relation arising from the three-nucleon interaction contain tensor-like terms correlating 
the spins and orientations of the nucleons. The contributions of these correlations to the 
response are discussed below. 

While these wave functions are not exact, they offer a rather precise characteriza- 
tion of the Euclidean response, as evidenced by comparisons with calculations using the 
correlated-hyperspherical-harmonics wave functions [flH] in A=3. These comparisons are 
presented in Sec. ||]. The Hamiltonian used in these studies is the Argonne model v' 8 
flTf N-N interaction plus the UIX three-nucleon interaction. This interaction reproduces 
many known properties of the alpha particle, including its binding energy and charge form 
factor. 

Calculation of the Euclidean response is a straightforward extension of the ground- 
state techniques employed in Green's function Monte Carlo. We wish to calculate matrix 
elements of the following type: 

(0\O 2 exp[-(H - E )r] O 1 \0) 



M(r) 



(0\exp[-(H-E )r)\0) 



(14) 
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For a ground-state calculation of the energy (Oi=l, O2—H) the matrix element is eval- 
uated by a Monte Carlo sampling of the coordinate-space paths. The denominator is 
exactly one for an exact ground-state wave function, otherwise there is a correction for 
finite t. For each path a complete set of 2 A A\/[N\ Z\) amplitudes is kept corresponding 
to all possible spin-isospin components of the ground state wave function. Since the op- 
erators do not in general conserve isospin we cannot use the most compact isospin basis 
used in ground-state calculations. 

For a more general matrix element M we simply keep another complete set of am- 
plitudes for each operator Ox, each set of amplitudes corresponding to the full operator 
acting on the ground state. The paths are sampled precisely as in the ground-state cal- 



culation J49|, and hence unaffected by the operators Oi, Oi- This allows us to calculate 
the response to a variety of operators (charge, current, different momenta, etc.) simulta- 
neously. 

We have found it computationally advantageous to calculate the response simultane- 
ously for several different directions of momentum transfer. A randomly picked set of 
three orthogonal axes are chosen, with q directions along both the positive and negative 
directions of each axis. This method yields much lower statistical errors in calculat- 
ing the response, and along with the more efficient methods for sampling path integrals 



recently applied to ground-state calculations [|{|, allows for much more precise results 
than obtained previously. It is also possible to calculate the response at several different 
momentum transfers simultaneously. 

It is certainly possible to extract more detailed information from the Euclidean re- 
sponse. Most efforts in this direction proceed under Maximum Entropy techniques em- 



ploying Bayesian statistics | 50fl . These techniques make use of the correlated error esti- 
mates in R(t) for different r. Given the enhanced precision of the present calculations we 
are exploring these possibilities. These considerations are beyond the scope of the present 
investigations, though, where we are primarily concerned with the total strength in the 
longitudinal and transverse channels. 



5 Electromagnetic current operator 

The model for the nuclear electromagnetic current adopted in the present study is briefly 
reviewed in this section for completeness, for a more complete description see Ref. [|5l] . 
The charge and current operators consist of one- and two-body terms: 

p(q) = EA (1) (q) + E^?(q) > ( 15 ) 

i i<j 

j(q) = Ei^fo) + £$(<*) > ( 16 ) 

i i<j 
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where q is the momentum transfer. The one-body operators p\ and j| have the standard 
expressions obtained from a relativistic reduction of the covariant single- nucleon current, 
and are listed below for convenience. The charge operator is written as 



pl 1) (q) = pl 1 iiR(q) + pSSc(q) , (17) 



with 



" So(q) = ' TTovw ~T e "" " ^ (2 " f,) q ' (at x P,) e """ ' (19) 



where Q 2 = q 2 — uj 2 is the four-momentum transfer, and u is the energy transfer. The 
current operator is expressed as 

i(1)(q) = [p * ' eiq ' r ' ] + " ^ q x a * giqri ' (20) 

where [■ • • , • ■ ■]+ denotes the anticommutator. The following definitions have been intro- 
duced: 

e< = G' fiiP (Q 2 )^(l + r Z)i ) + G' fiin (Q 2 )^(l-r Z)i ) , (21) 

IH = G M , p (g 2 )i(l + r Zil ) + G M ,n(Q 2 )^(l-r 2ii ) , (22) 

and p, a, and r are the nucleon's momentum, Pauli spin and isospin operators, respec- 
tively. The two terms proportional to 1/ m 2 in pj^c are the well known Darwin- Foldy and 
spin-orbit relativistic corrections |52], ^3], respectively. 

The calculations of the response functions discussed in the previous section have been 
carried out using the dipole parameterization of the nucleon form factors 

GeAQ 2 ) = Gd(Q 2 ) , (23) 

G E ,n{Q 2 ) = -/^rff^ , <"> 

GmAQ 2 ) = V P G D (Q 2 ) , (25) 
G M ,n(Q 2 ) = VnG D (Q 2 ) , (26) 

where ^ 

Gd(Q 2 ) = (1 + g2/A2)2 , (27) 

with A = 0.834 GeV/c, and where p p and p n are the proton (p p = 2.793 n.m.) and 
neutron (p n = —1.913 n.m.) magnetic moments, respectively. It is worth emphasizing 
that the available semi-empirical parameterizations of the proton electric and magnetic, 
and neutron magnetic form factors do not differ significantly — less than a couple of % 
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- in the low momentum transfer range of interest here, Q 2 < 0.4 (GeV/c) 2 , and that 
uncertainties in the neutron electric form factor have a negligible impact on the present 
results. Finally, we should note that in the actual calculations of the Euclidean responses 
the value of the four-momentum transfer occurring in the nucleon form factors (as well as 
in the electromagnetic iVA transition form factor, see below) is kept fixed at the quasi- 
elastic peak, as already mentioned in Sec. [4[ 

The most important features of the two-body parts of the electromagnetic current 



operator are summarized below. The reader is referred to Refs. |54|, |5T|| for a derivation 
and listing of their explicit expressions. 



5.1 Two-body current operators 



The two-body current operator consists of "model- independent" and "model-dependent" 
components, in the classification scheme of Riska p5| |. The model-independent terms are 
obtained [[5(J from the nucleon-nucleon interaction (the charge-independent part of the 
Argonne v is in the present study), and by construction satisfy current conservation with it. 
The leading operator is the isovector "7r-like " current derived from the isospin-dependent 
spin-spin (err) and tensor (tr) interactions. The latter also generate an isovector "p- 
like" current, while additional model-independent isoscalar and isovector currents arise 
from the isospin-independent and isospin-dependent central and momentum-dependent 
interactions. These currents are short-ranged and numerically far less important than the 
7r-like current. For the purpose of later discussions, we list below the explicit expression 
for the latter: 



4 2) (q; 



3)z 



e iq - ri /ps(r) Oi fa • r) + ^ f PS (r) a, fa • r) 



fa ■ V,)fa ■ Vj)(Vi - V i )^ PS (q; R, r) 



(28) 



where G\(Q 2 ) = Ge, p (Q 2 ) + GE,n(Q 2 ) is the isovector combination of the nucleon electric 
form factors, and R and r are the center-of-mass and relative positions of nucleons % and 
j, R = (r< + Tj)/2 and r = — Tj, respectively. The functions fp$ and gps are defined as 



dk 



Jk-r 



vps(k) , 



3ps(q; R,r) 



(2n) 3 S(q-k l -k j ] 



vps(ki) ~ vp S (kj 



p _ p 



(2tt) 3 (2vr) 3 

where vps{k) is obtained from the err and tr components of the interaction, 

vps(k) =v^(k)-2v tT (k) , 



(29) 
(30) 

(31) 
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with 



vaT ( k ) = tfl r 2 dr[ 3o (kr)-l]v^(r), (32) 



47T 



oc 



v tT (k) = — I r 2 drj 2 (kr)v tT (r) . (33) 

rv 







The factor jo(kr) — 1 in the expression for v aT (k) ensures that its volume integral van- 
ishes [061. 



In a one-boson-exchange (OBE) model, in which the isospin-dependent spin-spin and 
tensor interactions are due to 7r-meson (and p-meson) exchanges, the function vps(k) 
simply reduces to 

Vp8 ( k )^v«{k) = -£ T £f& 1 , (34) 
m~ k + m% 

where m n , f n , and f^ik) denote, respectively, the pion mass, ttNN coupling constant and 
form factor. In this limit, the functions fps and gps read: 

fps(r) - Mr) = ^7 ^(1 + m w r) (35) 

Arc {m 7T r) z 

iq-R ,.+1/2 e -L^(x)r 

g PS (q; R, r) - ^(q; R, r) = — — / dxe^ 1 -—— , (36) 

07T 7-1/2 L n {X) 

with 



L^x) = sJml + q^l-Ax^jA , (37) 

where for simplicity the irNN form factor has been set equal to one. The resulting current 
is then identical to that commonly used in the literature. 

The model-dependent currents are purely transverse and therefore cannot be directly 
linked to the underlying two-nucleon interaction. The present calculation includes the 
isoscalar pir'j and isovector uTc-f transition currents as well as the isovector current associ- 
ated with excitation of intermediate A-isobar resonances. The pir'j and uTcy couplings are 
known from the measured widths of the radiative decays p — > tcj [[jT] and uj — > [^H], . 
respectively, while their momentum-transfer dependence is modeled using vector-meson- 
dominance. Monopole form factors are introduced at the meson-baryon vertices with 
cutoff values of A 7r =3.8 fm" 1 and A p =A a; =6.3 fm" 1 at the ttNN, pNN and uNN ver- 
tices, respectively. 

Among the model-dependent currents, however, those associated with the A-isobar are 
the most important ones. In the present calculation, these currents are treated in the static 
A approximation rather than in the more accurate transition-correlation-operator scheme, 
developed in Ref. |6(J and applied to the calculation of the trinucleon form factors |n 



nd and pd radiative capture cross sections at low energies |54], £2), and S-factor of the 
proton weak capture on 3 He [ S3 1 . Again for later convenience, it is useful to list explicitly 
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the two-body A-excitation current used in the present work: 

, VNN^AN,ij 



Jg>< 



q;A) 



+ 



j,(q;A^iV)- 

m — rriA 

VAN^NN,ij . , AT a \ 

-j*(q;N -> A) 



m — m/\ 

where the iV ^ A electromagnetic current is modeled as 



(3f 



A) = ~G. 
2m 



7 iVA(g 2 )e iq - ri qxSiT, 



(39) 



and the expression for jj(q; A — > iV) is obtained from that for jj(q; iV — > A) by replacing 
the transition spin and isospin operators S and T with their hermitian conjugates. The 
electromagnetic form factor G 7 tva(Q 2 ) is parameterized as 



(40) 



where the N 



(1 + Q 2 /A%A,l) V 1 + Q 2 / A NA,2 ' 

A transition magnetic moment li* is taken here to be equal to 3 n.m., 
as obtained in an analysis of 7A 7 " data in the A- resonance region ] . This analysis also 
gives AjvA,i=0.84 GeV/c and AjvA,2=l-2 GeV/c. It is important to point out, however, 
that the quark-model value for /!*, jx* = (3a/2/ 5)/x]^ = 3.993 n.m. (/^ is the nucleon 
isovector magnetic moment), is often used in the literature. This value is significantly 
larger than that adopted above. Finally, the transition interaction vjy at— »AJV,y is given by 



VNN—>AN,ij 



K II (r)S i . ( r i + ^ n (r)4]T i .T : 



3 ' 



(41) 



and v an^n n ,ij is the hermitian conjugate of the expression above. The is the tensor 
operator where the Pauli spin has been replaced by the transition spin S«, and the 
functions v aTll (r) and v trll (r) are defined as 



v aTll (r) 



v tTa (r) 



f n f*m n e 
4/T 3 x 

4tt 3 



C{x) , 



X 



X 



(42) 
(43) 



where x = m^r, /* = (6y/2/5)f n is the quark-model value for the 7riVA coupling constant 
(adopted in the present work), and the cutoff function C(x) = 1 — e~ Xx2 , with A=4.09. 

Standard manipulations of the product of spin and isospin transition operators pD[ 
lead to the following expression for the A-excitation current: 



q;A) 



+ 



.G 7 tva(<3 2 
9m 

» ^ J , 



4r ; 



2>J 



U{r)<Tj + ?A(r)f(ffj • r) 



fA{r){oi x ctj) +^ A (r)(cTi x r)(crj 



x q 



(44) 
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where 



yarll /y\ ytrll/y,. 



f A (r) = ^ ^ , (45) 

m — niA 



v tTlI (r) 



g A (r) = 3 ±L . (46) 

m — mA 

The expression above reduces to that commonly used in the literature, if the quark-model 
values for the 7riVA and 7iVA coupling constants are adopted. 



5.2 Two-body charge operators 

While the main parts of the two-body currents are linked to the form of the two-nucleon 
interaction through the continuity equation, the most important two-body charge opera- 
tors are model-dependent, and should be considered as relativistic corrections. Indeed, a 
consistent calculation of two-body charge effects in nuclei would require the inclusion of 
relativistic effects in both the interaction models and nuclear wave functions. Such a pro- 
gram is yet to be carried out for systems with A > 3. There are nevertheless rather clear 
indications for the relevance of two-body charge operators from the failure of the impulse 
approximation (IA) in predicting the deuteron tensor polarization observable ||65|| , and 
charge form factors of the three- and four-nucleon systems pi], BH] . The model commonly 



used [67] includes the 7T-, p-, and cu-meson exchange charge operators with both isoscalar 
and isovector components, as well as the (isoscalar) pn^y and (isovector) loti^ charge tran- 
sition couplings, in addition to the single-nucleon Darwin- Foldy and spin-orbit relativistic 
corrections. The ir- and p-meson exchange charge operators are constructed from the 
isospin-dependent spin-spin and tensor components of the two-nucleon interaction (again, 
the Argonne vi& model), using the same prescription adopted for the corresponding cur- 
rent operators. Explicit expressions for these operators can be found in Ref. Here, 
we only emphasize that for Q < 1 GeV/c the contribution due to the 7r-exchange charge 
operator is typically an order of magnitude larger than that of any of the remaining 
two-body mechanisms and one-body relativistic corrections. 



6 Model studies 

The Euclidean response is an excellent tool to test our understanding of inclusive quasi- 
elastic scattering, since it incorporates an exact treatment of the states in the continuum. 
The Euclidean response does have the disadvantage, however, that it corresponds to a 
weighted integral over the energy loss ui and, as a consequence, the interpretation of 
potential differences between calculated results and experimental data is not so straight- 
forward. 
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Figure 5: Top: model response (solid line) changed by various mod- 
ifications (see text). Bottom: corresponding ratio of modified and 
unmodified Euclidean responses. 



In order to develop a better feeling for the properties of the Euclidean response, in 
this section we discuss a simple-minded model calculation. We use a parameterized cross 
section — de facto a fit to the longitudinal Rl{u) at one momentum transfer — and 
study the change in the Euclidean response upon various changes of the cross section as 
a function of energy loss. 
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In the top panel of Fig. [5] we show the model quasi-elastic peak as a function of energy 
loss (solid curve) and a selection of modifications. The changes have in general been 
made by adding a gaussian with arbitrary amplitude and selected position in energy loss. 
Figure |^ shows the quasi-elastic peak a) with a gaussian placed at very large energy loss 
(400 MeV), b) a gaussian placed on the high-energy loss tail of the quasi-elastic peak (200 
MeV) and c) a gaussian placed on the low-o; side of the peak (100 MeV). It also displays 
a curve where d) the width of the quasi-elastic response has been decreased by 20% (with 
the overall amplitude adjusted to conserve the area) and one for the case where e) the 
quasi-elastic peak is shifted by 10 MeV. 

The lower panel of Fig. [5] shows the resulting changes in terms of the ratio of modified 
to original Euclidean responses. The value at r=0 reflects the (arbitrary) integral over uo 
of the added modification. Several features are noteworthy: 

• The Euclidean response at finite r very quickly suppresses the contribution from 
large energy loss. The dash-dot curve shows that already at r > 0.01 the contribu- 
tion from the large peak added at uo=A00 MeV is suppressed. For the experimental 
transverse response function Rt{uo) this implies that the contribution from pion pro- 
duction in the A peak (which is not included in the theory we are going to compare 
to) is only affecting the results for very small r. We will therefore ignore this region. 

• The region of the quasi-elastic cross section at low uo comes in very prominently at 
the larger values of r, as indicated by the curve labeled "peak at 100 MeV" . 

• A shift of the quasi-elastic peak to larger uo leads to an Euclidean response that 
quickly falls with increasing r, reaching saturation by the time r gets to values 
approaching 0.05. 



7 Results 

We have used the longitudinal (L) and transverse (T) experimental response functions 
of Figs. [I] and to compute the corresponding experimental Euclidean responses shown 
in Figs. |6] and ^. The nucleon electromagnetic form factors are divided out using the 
parameterizations of Holder et al. | 68[| . In order to not include too much of the tail of the 



A-resonance, the integration has been performed up to the energy loss uo where the T- 
response starts to increase significantly with uo (the corresponding value of uo is indicated 
in Figs. [1] and ^ by a +). Since for the T-Euclidean response at very small r the tail of the 
A-peak nevertheless plays a role, the experimental response in this region is indicated by 
a dashed line only, and should not be compared to the theoretical calculations discussed 
below. 

The statistical errors of the experimental Euclidean response are obtained via the 
usual error propagation when integrating. The additional overall systematic uncertainty, 
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Figure 6: Longitudinal (upper half of figure) and transverse Eu- 
clidean response of 3 He for momentum transfers 300-600 MeV/c. 
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Figure 7: Longitudinal (upper half of figure) and transverse Eu- 
clidean response of 4 He for momentum transfers 300-600 MeV/c. 
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estimated from the L/T-separation, amounts to typically 3% for 3 He for both the L- and 
T-cross sections. For 4 He a similar uncertainty in the scale factor applies at the lower 
momentum transfers and for both L and T; at 600 MeV/c the uncertainty in Rl{u) 
increases to 6%. These scale errors, which then apply also to the corresponding Euclidean 
responses, have not been included in the error bars shown in Figs. |^ and |7|. 

In Figs. H and [7| we also show the calculated Euclidean responses, obtained in IA 
and with inclusion, in addition, of the contributions associated with the two-body charge 
and current operators, discussed in Sec. ^|. It is immediately apparent that two-body 
contributions reduce by a small amount the L-responses, while increasing the T-responses 
very substantially at all momentum transfers. The enhancement in the T-channel occurs 
already at low u, as is seen from the Euclidean response at large r. Two-body effects thus 
are important over the entire quasi-elastic peak, and not only — as was often expected 
- in the "dip- region" on the large-a; side of the quasi-elastic peak. These conclusions are 
in agreement with those of an earlier study [^T|, as well as with those inferred from the 



super-scaling analysis of Ref. 0] for nuclei with mass number A— 12-56. 



When considering the effect of two-body currents as a function of momentum transfer 
- in particular, when studying Fig. |7| — one notes that at low q the effect of two-body 
currents at large to (low r) is bigger than at low u (large r). At large q, this situation 
becomes the reverse. Figures |6| and [7] also show that theory explains well the rapid 
increase of two-body contributions between 3 He and 4 He. In contrast to most published 
calculations (for a discussion see Sec. |l]), the present calculation does give the sizeable 
two-body contribution required by the data. 

Figure [7] shows that the T-Euclidean response at low q rises very rapidly towards very 
small r, reaching almost twice the IA value at r=0, thus suggesting that part of the two- 
body strength is located at very large u, basically under the A-peak (compare to Fig. |5|). 
It also implies that this strength is very spread out in uj, and presumably best discussed 
in terms of the sum-rule (see Sec. ||). 

At lower values of momentum transfer, the calculated 4 He T-Euclidean response is a 
bit high at large r, implying that the corresponding calculated cross section would be 
somewhat too high at low u. As emphasized by the sensitivity studies in Sec. |], the low 
uj region gets great weight for large r, so a small increase in the absolute value of o{uS) 
leads to a large increase in E(t). 

Overall, the agreement between theory and experiment for 4 He, the nucleus which 
allows us best to study the relative role of one- and two-body contributions, is excellent 
for the L-response, thus implying that an accurate treatment of the nuclear spectrum has 
been achieved, since two-body operators give small corrections in the L-channel. For the 
4 He T-Euclidean response, the large two-body effects predicted by theory are confirmed 
by experiment, although the associated contributions are a bit too large in the g-range 
400-500 MeV/c. 
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8 Longitudinal and transverse sum rules 



Sum rules provide a powerful tool for studying integral properties of the response of the 
nuclear many-body system to an external electromagnetic probe. Of particular interest are 
those for the longitudinal and transverse response functions at constant three-momentum 
transfers, since they can be expressed as ground-state expectation values of the charge 
and current operators and, therefore, do not require any knowledge of the complicated 
structure of the nuclear excitation spectrum. Direct comparison between the theoreti- 
cally calculated and experimentally extracted sum rules cannot be made unambiguously, 
however, for two reasons. Firstly, the experimental determination of the longitudinal and 
transverse sum rules requires measuring the associated response functions in the whole 
energy-transfer range, from threshold up to oo. Inclusive electron scattering experiments 
only allow access to the space-like region of the four-momentum transfer (u < q). While 
the response in the time-like region (u > q) could in principle be measured via e + e~ 
annihilation, no such experiments have been carried out to date, to the best of our knowl- 
edge. Therefore, for a meaningful comparison between theory and experiment, one needs 
to estimate the strength outside the region covered by experiment. In the past, this has 
been accomplished, in the case of the longitudinal response, either by extrapolating the 
data or by parameterizing the high-energy tail and using energy-weighted sum rules 
to constrain it [[TJ], |7T[. For the A=2-4 nuclei, the unobserved strength amounts to 5-10% 
at the most for three-momentum transfers in the range q < 1 GeV/c ||71|| , and both pro- 
cedures lead to similar results. Indeed, the calculated (non-energy-weighted) longitudinal 
sum rule — also known as the Coulomb sum rule — appears to be well satisfied by the 
data |7g,[7g. 



The second reason that makes the direct comparison between theoretical and "experi- 
mental" sum rules difficult lies in the inherent inadequacy of the present theoretical model 
for the nuclear electromagnetic current, in particular its lack of explicit pion production 
mechanisms. The latter mostly affect the transverse response and make its A-peak region 
outside the boundary of applicability of the present theory. The charge and current oper- 
ators discussed in Sec. |j, however, should provide a realistic and quantitative description 
of both longitudinal and transverse response functions in the quasi-elastic peak region, 
where nucleon and (virtual) pion degrees of freedom are expected to be dominant. In 
light nuclei and at the momentum transfer values of interest here, the quasi-elastic and 
A-production peaks are well separated, and it is therefore reasonable to study sum rules 
of the transverse response. 

While non-energy- and energy-weighted longitudinal sum rules have been extensively 
studied in the past (see Refs. |ETL 1731 for a review), the number of studies dealing with sum 



rules of the transverse response is much more limited [74]. The present section focuses 
on the latter, in particular on the enhancement of transverse strength due to many-body 
components of the electromagnetic current, within the limitations discussed above. It 
also addresses, within the sum-rule context, the issue of the enhancement in the ratio 
of transverse to longitudinal strength, observed in the quasi-elastic response functions of 
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nuclei. Finally, it attempts to provide a semi-quantitative explanation for the observed 
systematics in the excess of transverse strength, both as function of mass number and 
momentum transfer. All the calculations are based on the AV18/UIX Hamiltonian model, 
and use correlated-hyperspherical-harmonics (variational Monte Carlo) wave functions for 
A=3-A (A=6) nuclei. 

The (non-energy-weighted) sum rules are defined as 
S a (q) = C a du;S a (q,uj) 

Ju} th 

= C Q [(0|Ot(q)O Q (q)|0)-|(0|O Q (q)|0)| 2 ] , (47) 

where S a (q, uj) is the point-nucleon longitudinal (a=L) or transverse (a=T) response func- 
tion, O a (q) is either the charge p(q) or current j(q) operator divided by the square of the 
proton form factor \G P E (Q 2 )\ 2 (again, Q 2 is evaluated at the energy transfer corresponding 
to the quasi-elastic peak), |0) denotes the ground state, and the elastic contribution to 
the sum rule has been removed. An average over the nuclear spin orientations is tacitly 
implied in the evaluation of the expectation values. The constant C a , for a=L or T, is 
given by 

Cl = \ , (48) 
2 m 2 1 

where Z (N) and fi p (fi n ) are the proton (neutron) number and magnetic moment, re- 
spectively. It has been introduced in Eq.(f|7]) so that, in the limit q — > oo and under the 
approximation that the nuclear charge and current operators originate, respectively, from 
the charge and spin-magnetization of the individual nucleons only, S a (q — > oo) = 1. Note 
that the Euclidean response functions calculated in Sec. [I] and the sum rules defined here 
are related via 

S a (q) = C a E a (q,r = 0) . (50) 

The expectation values in Eq. (E7p are calculated with Monte Carlo methods, without 
any approximations. 

The calculated sum rules for 3 He, 4 He, and 6 Li are listed in Tables |I] and 0. The 
longitudinal sum rule S^q) is relatively un-influenced by two-body charge operators, in 
agreement with the results of an earlier study ||72|| . The transverse sum rule S'r(g) is 
substantially increased by two-body current contributions. The resulting enhancement 
has two interesting features: i) it increases, for fixed q, in going from A=(2 to) 3 to 4, and 
decreases from A=4 to 6; ii) it decreases, for fixed A, as q increases. Both these features 
are summarized in Figs. || and |9|, in which the ratios Sx(q) / S^q), obtained by including 
one-body only and both one- and two-body contributions, are plotted as function of A 
for fixed q and as function of q for fixed A. The former figure is reminiscent of Fig. |], in 
which the ratio of transverse to longitudinal strength in the quasi-elastic region is obtained 
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3 He 4 He 6 Li 


q(MeV/c) 1 1+2 


1 1+2 


1 1+2 


300 0.787 0.763 
400 0.921 0.875 
500 0.964 0.901 
600 0.982 0.908 
700 0.994 0.914 


0.670 0.649 
0.859 0.815 
0.941 0.881 
0.973 0.910 
0.994 0.942 


0.977 0.933 
0.995 0.932 
0.990 0.921 
0.990 0.924 
0.994 0.938 



Table 1: The longitudinal sum rule obtained with one-body only and both one- and 
two-body charge operators. 



3 He 4 He 6 Li 


q(MeV/c) 1 1+2 


1 1+2 


1 1+2 


300 0.929 1.31 
400 0.987 1.30 
500 1.01 1.28 
600 1.01 1.25 
700 1.01 1.23 


0.893 1.67 
0.970 1.62 

1.00 1.55 

1.01 1.49 
1.01 1.44 


0.912 1.57 
0.974 1.52 
0.999 1.46 
1.01 1.41 
1.011 1.37 



Table 2: The transverse sum rule obtained with one-body only and both one- and 
two-body current operators. 



from the measured response functions. Obviously, the truncated integrals in Fig. |] do not 
include the strength at high uj. 

The purpose of the present section is to offer an explanation of the features mentioned 
above. To this end, three points are worth emphasizing. Firstly, among the two-body 
current contributions, the most important are those associated with the PS (pion-like) 
and A-excitation currents. This fact has been explicitly verified by direct calculation, 
as shown in Table |3| for 4 He, as an example. Note that the results in the 2nd and 3rd 
columns are slightly different from those reported above in Table since they are based 
on a random walk consisting only of 1,000 configurations, much shorter than that used in 
the calculations of Table §. These calculations, though, are based upon the same random 
walk and therefore allow a better determination of the individual contributions. 

Secondly, consider expanding the current into one- and two-body components ji and 

3lmi 

J = Y^^l + Y^^lm ■ (51) 

I l<m 

Then, ignoring the very small (and, with increasing q, rapidly vanishing) elastic contri- 
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Figure 8: The ratios St{q)/ Sl(q), obtained with one-body currents only and both one- 
and two-body currents, as function of mass number A. 



q(MeV/c) 


1 


1+2 


1+2-7T+2-A 


300 


0.915 


1.65 


1.58 


400 


0.980 


1.59 


1.50 


500 


1.01 


1.53 


1.44 


600 


1.01 


1.47 


1.38 


700 


1.01 


1.41 


1.33 



Table 3: The 4 He transverse sum rule. 



bution to St{q), one finds that the first term in Eq. fl4"7|) 

+ X! [0'J + jln)jlm + h.C.] + X Jlmjlm 
l<m l<m 

+ terms involving 3 or 4 different nucleons . (52) 
At large momentum transfers, one would expect terms involving 3 or 4 nucleons to be 
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Figure 9: The ratios St{q)/ Sl(q), obtained with one-body currents only and both one- 
and two-body currents, as function of momentum transfer q. 

small, particularly in light nuclei where Pauli correlations are unimportant. Dropping the 
last term corresponds to considering only incoherent scattering from pairs of nucleons. 

This simple expectation is indeed borne out by a direct calculation, the results of 
which are listed for 4 He in Table |j. Thus the transverse sum rule appears to be saturated 
by the one- and two-body terms in the expansion for j*j above. 



q(MeV/c) 


1 


1+2 


1+2-reduced 


300 


0.915 


1.65 


1.70 


400 


0.980 


1.59 


1.59 


500 


1.01 


1.53 


1.51 


600 


1.01 


1.47 


1.45 


700 


1.01 


1.41 


1.39 



Table 4: The 4 He transverse sum rule: effect of three- or four-nucleon terms. 



Thirdly, the transverse strength associated with two-body currents is almost entirely 
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due to pn pairs. To make this observation more precise, consider the "reduced" two-body 
current: 

Jim * Jim 

{P l P m + N l N m ) , (53) 

where Pi and iVj are the proton and neutron projection operators for particle I. Thus 
the "reduced" two-body current only acts on pp or nn pairs, and the transverse sum rule 
calculated with it should be given almost entirely by the one-body part of j. This fact 
is again confirmed by direct calculation, as it is evident from Table [5[ That pn pairs 



q(MeV/c) 


1 


1+2 


1+2-pp or nn only 


300 


0.915 


1.65 


0.919 


400 


0.980 


1.59 


0.987 


500 


1.01 


1.53 


1.02 


600 


1.01 


1.47 


1.03 


700 


1.01 


1.41 


1.03 



Table 5: The 4 He transverse sum rule: contribution of pp and nn pairs. 



are responsible for the strength due to two-body currents can also be understood by 
the following considerations. The pion-like and A-excitation currents have the isospin 
structure (see Sec. |5|), again in a schematic notation, 

jim{n) = {n x T m ) z Oi m (n) , (54) 

ji m (A) = Ti jZ Oi m (A, a) + r mi2 mi (A, a) + (77 x r m ) 2 0; m (A, b) , (55) 
while the leading part of the one-body current is given by 

ji = t^O^IV) , (56) 

where Oi(IV) denotes the isovector part of j t . Now the term jj m ji m (with j im including 
pion-like and A-excitation currents) will produce, as far as isospin is concerned, terms 
like 

(ti x T m )l = 2(1- T hz T m , z ) , (57) 
(Tl,z Or T m>z )(Ti x r m ) z = ±i (n ■ r m - T l)Z T m)Z ) , (58) 

_ T\m + Ti-T m 

<l, Z 1 m,z — g j W y / 



where T lm is the isotensor term T/ m = 3r; z r m 2 — 77 ■ r m . 

In addition, there will also be isospin-independent terms, of the type 

OL(A,a)O im (A,a) . (60) 



However, it is important to note that the operators in Eqs. (57) and (pq) vanish when 



acting on pp or nn pairs. It should also be noted that the isotensor term in Eq. 

28 



vanishes in T=0 and T=l/2 ground states, namely in 3 He, 4 He and 6 Li. A similar analysis 
can be carried out for the interference terms between one- and two-body currents, 

J2(jj + jl)jlm + h.C. , (61) 
l<m 

for which one obtains isospin-independent, and type ( |5g|) or type (|59|) operators. In any 
case, the direct calculation indicates (see Table |5]) that pp and nn pairs do not contribute 
appreciably to Sriq)- 

On the basis of the above observations and ignoring the convection term in the one- 
body ji, one concludes that the excess transverse strength, defined as 

AS T (q) = S T (q) - S T (q; 1-body) , (62) 

must be proportional to 

POO 

ASr(q) = / dxtr[F(x; q) p(x;pn)] , (63) 
Jo 

where F(x; q) is a complicated matrix in the spin-space of the two nucleons depending 
upon the current operators alone, and the A-dependence is included in the pn elements of 
the two-nucleon density matrix p 2 (x;pn, sj, s m , s' t , s' m ). Here si, s m , etc., are spin projec- 
tions (up or down) of particles /, m, etc. In fact, one can express these densities in terms 
of total spin-isospin S, T=0,1 or 1,0 for pair lm. The crucial point is that, in nuclei, these 
pn densities scale, see |75| , namely 

p 2 (^; T = 0; A) = R A p(x; T = 0; deuteron) , (64) 

p 2 (^; T = 1; A) = R' A p(x; T = 1; S quasi— bound) . (65) 

The scaling factors Ra and R' A have been calculated in Ref. |75| , with R' A ~ R A and 
i?A=2.0, 4.7, 6.3, 18.8 for 3 He, 4 He, 6 Li, and 16 0, respectively, and so one would expect 
ASx(q) to scale with 

oc , (66) 

where the factor in the denominator on the r.h.s. is from the normalization adopted for 
S T (q). 

The calculated values for the excess strength ASV(?) are listed in Table []. On the 
basis of the scaling law above one would deduce 

|g - 0,67 1. 0, 9 4, (68, 
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q(MeV/c) 


3 He 


4 He 


6 Li 


300 


0.38 


0.78 


0.66 


500 


0.27 


0.55 


0.46 


700 


0.21 


0.43 


0.36 



Table 6: The excess strength ASV(<?) calculated in 3 He, 4 He, and 6 Li. 



and these values are reasonably close to those of Table ^| They are also close to those that 
can be inferred from data, see Fig. Finally, in Fig. |IU] the integrands in Eq. ( |B"3"D are 
displayed for 3 He, 4 He, and 6 Li, properly scaled according to the factor in Eq. ([661) . Note 
that also shown are the contributions due to pn pairs in T=0 states only. The behavior 
of the integrands, as illustrated in Fig. |H], is to be expected, since it is a consequence 
of the "scaling" behavior more generally observed for the calculated T, 5=0,1 and 1,0 
pair-distribution functions in nuclei |75| . 




Figure 10: The integrands in Eq. (|63"D for 3 He, 4 He, and 6 Li at momentum transfers of 
300 and 700 MeV/c, scaled according to the factor in Eq. (|66"D. Also shown are the 
contributions due pn pairs in T=0 states only. 

Note that the dominant contributions to the excess strength occur for pair separations 
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around or slightly less than one fermi. One would naturally associate this strength with 
a significant contribution to two-nucleon final states of relatively large relative momenta. 
Detailed microscopic calculations in A=4 with full final-state interactions and two-body 
currents will be necessary to make precise predictions. 

It should be emphasized that the scaling law for ASV(<?) can be used to estimate the 
excess transverse strength in nuclei, once the factors Ra are known. In nuclear matter, for 
example, the authors of Ref. [TBI obtain Roo = 1.59, the latter being defined as Ra/A in 



the limit A — > oo, and therefore one would expect a very substantial enhancement of the 
transverse sum rule (and, consequently, transverse response function) due to two-body 
currents, namely ASr{q] oo)/ASV(<?; 4 He) ~ 1.35. 

The nuclear matter transverse response has been calculated in Ref. |77]], by using 
correlated-basis-function perturbation theory and including, in addition to the single- 
nucleon spin- magnetization current, the pion-like and A-excitation two-body currents. 
Explicit integration of the response functions |77j, [78] indicates that the transverse strength 
is increased by the two-body contributions by roughly 15% in the momentum transfer 
range 300-700 MeV/c. This enhancement is significantly smaller than that inferred from 
the scaling law above. The underestimate is presumably due to the inherent limitations 
in the calculations carried out so far, which only retained one-particle-one-hole (Ip-lh) 
intermediate states and estimated the contribution of two-particle-two-hole states by 
folding the Ip-lh response with a width derived from the imaginary part of the optical 
potential. It should be possible to calculate the transverse sum rule by direct evaluation 



of the ground state expectation value. Work along these line is in progress |78| 



The longitudinal and transverse sum rules in matter can be estimated in a Fermi gas 
model in a similar simple manner. As in calculations of 4 He (see Sec. ||), this is useful 
to help understand the role of initial-state correlations in the transverse response of the 
nucleus. We again ignore contributions of 3 and 4 nucleon terms as in Eq. [52[ In matter 
this approximation should be valid at high momentum transfer q, but becomes more 
questionable as q is decreased. A significant enhancement of the transverse sum rule is 
expected due to the short-range part of the two-body currents — these necessarily involve 
large momenta between the pair of nucleons, thus broadening the range of validity of this 
simple approximation. 

The Fermi-gas sum rules are decomposed into parts depending only on the single 
nucleon currents and the remaining terms which also involve two-nucleon currents: 

S a (q) = S a (q; 1-body) + AS a (q) , (69) 

where a=L, T. The simplest approximation to the response due to one-body currents is to 
assume incoherent scattering from isolated nucleons. This yields S a (q) = 1, neglecting the 
neutron charge and convection current contributions to the longitudinal and transverse 
response functions, respectively. 

The additional contributions to the sum rules involving two-nucleon currents can be 
written as a sum of interference terms between one- and two-nucleon operators and the 
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square of the two-nucleon operators, in the same approximation adopted above. The 
short-range nature of the two-nucleon operators implies that incoherent scattering from 
pairs of nucleons should be dominant: 

AS a (q) — (A — 1)[2 ($|[C^(q)C^ m (q) + h.c] + 0^ m (q)O aM (q)|$}] , (70) 

where the factor 2 in the interference term arises because the pair term can connect to 
either of the two single-nucleon operators. 

To make a simple estimate of the contribution of the two-nucleon terms of the current, 
we consider the two-nucleon density matrix p2(ri m ) XtX ', which depends only upon the 
separation between the pair of nucleons and upon their initial and final spins and isospins 
X and x' ■ m the Fermi gas approximation, p 2 is diagonal in the spins and isospins, and 
the spatial dependence is given by simple Slater functions. We then obtain: 

ASM) =PE j j 1 - E P 2 (^) Xl x'(xl[2 0l(q) + Ol M (q)]0, m (q)ix / ) , (71) 

x,x' 

where momentum-dependent pieces in the current have again been dropped. The excess 
contributions involving two-nucleon currents are given in Table |7| for the Fermi gas model. 
The longitudinal contributions are positive but small, ranging up to ~ 0.02. The trans- 
verse contributions, ranging from ~ 0.06 at 700 MeV/c to ~ 0.11 at 300 MeV/c, are 
significant. 

Finally, as far as the g-dependence is concerned, from the explicit expressions of the 
current operators in Sec. |5] it is evident that the excess transverse strength should behave 
as 

AS T {q)~{oc + (3q + 1 q 2 )/q 2 , (72) 

where the g 2 -factor in the denominator is due to the normalization adopted for Sriq), 
and so will approach a constant in the limit of large q. 



q (MeV/c) 


ASl 


AS T 


300 


0.004 


0.114 


400 


0.007 


0.081 


500 


0.011 


0.066 


600 


0.017 


0.060 


700 


0.024 


0.056 



Table 7: Excess-strength contributions ASl and ASt to the Fermi gas sum rules from 
terms involving two-nucleon currents. 
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9 Model studies with simplified interactions, wave 
functions, and currents 



As a guide to better understanding these results and comparing with other calculations, 
it is useful to contrast the complete calculations described above with various truncations 
of the initial ground state wave function, the current operators, and the Hamiltonian. Of 
course only the complete calculations can be meaningfully compared to the data, as they 
include both initial state wave functions and current operators which are consistent with 
the Hamiltonian used to determine the r- (or energy-) dependence of the response. 

The transverse channel is most interesting in this regard, as it shows a large enhance- 
ment from the two-nucleon current operators. Results for 4 He at 400 MeV/c with various 
truncations are shown in Fig. |TT|. The truncations include full (FW) and simple (SW) 
wave functions, full (FC) and impulse (IC) currents, and full (FI) and simple (SI) inter- 
actions. The simple wave functions and interactions are described below. The differences 
in the longitudinal channel are much less dramatic. 

The full ground-state variational wave function is described above (Eq. |T^ ), it includes 
strong tensor correlations from the pair correlation operators Fij and from the three- 
nucleon correlation [T^. In order to better determine the origin of the enhancement 
arising from the two-body currents, we have also considered a simplified ground-state 
wave function (SW) where the tensor correlations ii*(r) and u tT (r) (Eq. [13])) and the U^-k 
correlations arising from the two-pion-exchange three-nucleon interaction have been set 
to zero. 

Similarly it is interesting to compare the effect of different Hamiltonians describing 
the final-state interactions. In the Euclidean response this corresponds to using different 
Hamiltonians for the imaginary-time propagation of the system. The Hamiltonian used 
in the propagation does not directly affect the sum rules which depend only upon the 
initial state and the current operators. We have constructed a simplified V4 interaction 
(SI) where the tensor terms in the full Hamiltonian have been set to zero. This would, of 
course, yield a very under-bound alpha-particle ground-state. To compensate, we add a 
potential of two-pion exchange range to both the spin one channels: 

v' s=1 . T=0>1 (r) = v s= i;T=oAr) ~ 1-4 T 2 (r) , (73) 

where 



T 



1 — exp(— cr 2 ) 



e -^ t ^ 1 ( 74 ) 

/ir {fir) z J r 

is taken from the Argonne interaction models and is a function of two-pion exchange range. 
The constant 1.4 MeV-fm has been set to crudely reproduce the alpha-particle binding. 
This allows us to concentrate on the spin-dependence of the final-state interactions as 
opposed to drastically altering the spectra of the struck nucleus. In all cases the full 
currents (FC) are those obtained from the Argonne vis interaction (AV18), they have not 



33 



been reconstructed to be consistent with the Hamiltonian used for the initial- or final- 
states. The motivation here is to examine the various contributions to the full calculation. 



From the figure it is clear that a dramatic difference remains between full (FC) and 
impulse (IC) currents whatever model is chosen for the wave function and Hamiltonian. 
The t = (sum-rule) difference between full and impulse currents is largest for the full 
wave function, but even with a highly simplified wave function (SW) a large difference 
remains between the results with full currents (FC, SW) and impulse currents alone (IC, 
SW). 

On the basis of Fermi-gas calculations of matter, it had been believed that the large 
enhancement from two-body currents found in previous calculations of light nuclei [21] were 
due to the presence of strong tensor correlations in the ground-state wave function. While 
these correlations do make a significant contribution, even simplified wave functions show a 
substantial enhancement. In light nuclei, at least, this is a consequence of the complete set 
of final states automatically included in the sum rule and Euclidean response calculations. 
We have also considered more drastic simplifications for the ground-state wave function, 
including only central (f c ) correlations. Even in this case there is a dramatic enhancement 
of the response when two-nucleon currents are included. 

Of course the Hamiltonian used for final-state interactions cannot affect the sum rule at 
r=0, but it can change the energy- dependence of the response. Calculations with simple 
wave functions (SW) and the simplified interactions (SI) are shown as diamonds in the 
figure. With simplified wave functions, interactions and impulse currents (IC, SW, SI), the 
slope at r=0 is much more shallow corresponding to an energy-weighted sum rule much 
closer to k 2 / (2m) than in the full calculation — of course, this is to be expected, since 
tensor components, missing in the SI model, substantially enhance the energy-weighted 
sum rule. This same interaction also has a larger low-energy (large-r) response than the 
calculations made with the full current. This is undoubtedly related to the choice of 
modified Hamiltonian, the choice made here will be more attractive in p-waves than the 
full Hamiltonian, and these presumably dominate the low-energy transverse response in 
4 He. 

With the full currents (FC), there is much less dependence upon the choice of final- 
state interactions. Indeed, the calculations with the simple wave function and full currents 
(FC, SW) are nearly identical over the range of r considered. The low-energy p-wave 
continuum in the more attractive simplified Hamiltonian yields less overlap with the two- 
nucleon current operators, resulting in a very similar full response for the two different 
final-state interactions. The full calculation (FC, FW, FI) has a much larger contribution 
at higher energy, resulting in a steep initial fall-off with r. It also a has a somewhat 
smaller response at low energy (large r) than the full calculation. 

Finally, we have calculated the responses in A=3 using the correlated-hyperspherical- 



harmonics (CHH) wave functions obtained by the Pisa group [48[ for this same Hamilto- 
nian. Calculations of the longitudinal response of 3 He at various momentum transfers are 
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Figure 11: Euclidean transverse response for 4 He at 400 MeV/c with full (F) or IA (I) 
current (C), full or simplified (S) wave function (W), full or simplified interaction (I). 



compared in Fig. |12|. The differences between the variational and CHH wave functions 
are very small, as is apparent in the figure. This is perhaps not surprising, as the drastic 
truncations made in the comparisons of simple (SW) and full (FW) variational wave func- 
tions were themselves somewhat modest. Differences in the CHH and VMC transverse 
response calculations of 3 He are also quite small. 

These calculations demonstrate that the two-nucleon currents play a crucial role in the 
transverse response. Precise comparisons with experimental data also require calculations 
with accurate initial-state wave functions and final-state interactions. In such realistic 
calculations, the contributions of the two-nucleon currents are large both in the integrated 
response and in the low-u; regime. 
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Figure 12: Comparison of the 3 He Euclidean longitudinal response functions using varia- 
tional (open symbols) and CHH (solid symbols) wave functions. 



10 Conclusions 



We have determined the 3 He and 4 He longitudinal and transverse response functions in the 
momentum transfer range 300-700 MeV/c from an analysis of the (e,e') world data. The 
corresponding Euclidean response functions have been derived by direct Laplace trans- 
form of the experimental data, and have been found to be in satisfactory agreement with 
those calculated with Green's function Monte Carlo methods using realistic interactions 
and currents. Leading terms of the two-body charge and current operators are constructed 
consistently with the two-nucleon interaction included in the Hamiltonian (the Argonne 
Vis). A number of improvements in the algorithms employed for the Monte Carlo eval- 
uation of the relevant path integrals have allowed us to reduce, very significantly, the 
statistical errors in the Euclidean response calculations. 

Two-body charge operators reduce slightly the one-body longitudinal strength at large 
r (corresponding to the threshold region of Ri(q,uj)), while two-body currents increase 
very substantially, and particularly for 4 He, the one-body transverse strength over the 
whole r-range considered. Thus, in the quasi-elastic region, single-nucleon knock-out 
processes are dominant in the longitudinal channel, while both one- and two-body mech- 
anisms contribute with comparable magnitude in the transverse channel. These qualita- 
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tive conclusions are corroborated by the scaling analysis of the data described in Sec. [| 
the longitudinal and transverse scaling functions fi and fx, which would be expected 
to overlap if one-body processes alone were to be at play, display in fact drastically dif- 
ferent trends (see Fig. §). The enhancement in the ratio of transverse to longitudinal 
quasi-elastic strength can be quantified by considering integrals of fi and fx (of course, 
over the quasi-elastic peak region alone), as done in Fig. |j. Experimentally, this T/L 
ratio is found to increase very significantly from A=3 to 4, to decrease only moderately 
from A=A to 12, and to remain rather flat as A— 12-56. Of course, the interpretation of 
the integral of fa as reflecting exclusively quasi-elastic strength is not entirely correct, 
particularly since, as the momentum transfer becomes larger and larger, the quasi-elastic 
and A peaks tend to merge together: strength from the pion-production region will then 
necessarily spill over into the quasi-elastic region, contaminating fr- Nevertheless, the 
amount of "spurious" (non quasi-elastic) strength contamination should not be too large, 
at least for light nuclei, for which the quasi-elastic and A peaks remain well separated at 
all momentum transfers considered here. The observed enhancement of the T/L ratio in 
3 He and 4 He is well reproduced by theory, since the Euclidean response functions derived 
from data are close to those obtained in the calculations, over the whole r-range. 

The T/L ratio has also been studied in the A=3, 4, and 6 nuclei via sum-rule tech- 
niques. Even within the limitations that such an approach necessarily entails (see Sec. |8|), 
there are rather clear indications that the present theory is able to predict its observed 
dependence upon both mass number and momentum transfer, see Figs. |8| and ^. The sum 
rule study in Sec. |8| has also allowed us to establish quantitatively that the excess trans- 
verse strength associated with two-body currents is almost entirely due to pn pairs. This 
fact has then led to the scaling law for the excess transverse strength ASt{q) — Ra/A, 
which derives from the universal scaling behavior obtained for the calculated pn-pair 



distribution functions in nuclei 75 



Finally, the role of tensor interactions and correlations has been investigated via model 
studies of the 4 He Euclidean transverse response function, using simplified interactions, 
currents, and wave functions. In contrast to earlier speculations [ET] that the large en- 
hancement from two-body currents was due to the presence of strong tensor correlations in 
the ground state, it is now clear that this enhancement arises from the concerted interplay 
of tensor interactions and correlations in both ground and scattering states. A successful 
prediction of the longitudinal and transverse response functions in the quasi-elastic region 
demands an accurate description of nuclear dynamics, based on realistic interactions and 
currents. 
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